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We investigate the influence of electron-phonon interactions on the dynamical properties of a 
quantum-dot-cavity QED system. We show that non-Markovian effects in the phonon reservoir lead 
to strong changes in the dynamics, arising from photon-assisted dephasing processes, not present in 
Markovian treatments. A pronounced consequence is the emergence of a phonon induced spectral 
asymmetry when detuning the cavity from the quantum-dot resonance. The asymmetry can only be 
explained when considering the polaritonic quasi-particle nature of the quantum-dot-cavity system. 
Furthermore, a temperature induced reduction of the light-matter coupling strength is found to be 
relevant in interpreting experimental data, especially in the strong coupling regime. 
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The emergent field of quantum information technol- 
ogy [l[ has spurred major research activities on control- 
ling the fundamental interaction between a semiconduc- 
tor quantum-dot (QD) and a cavity. Solid-state cavity 
QED (cQED) systems are inherently coupled to the envi- 
ronment, since the emitter is embedded in a solid. This 
is in contrast to atomic cQED where the atom can be 
effectively isolated and only few discrete energy levels 
are sufficient in the description. Remarkably dephasing 
from solid-state environments cannot simply be seen as 
a nuisance, but can in fact lead to enhanced coupling of 
QDs to a detuned cavity mode of importance for efficient 
single-photon sources [l-Il] and nanolasers 0. Modeling 
the continuum of reservoir modes of solid-state systems 
constitutes a considerable challenge. The coupling of the 
QD-cavity system to its solid-state environment has al- 
most exclusively been described using Markovian theories 
0,0], neglecting memory effects of the reservoirs. While 
the Markovian approximation is well justified for some 
reservoirs, this is not in general true for the reservoir 
consisting of quantized lattice vibrations. Such phonon 
reservoirs dephase the QD-cavity system, whereby the 
entanglement between light and matter in general is de- 
stroyed. Notably the first experimental demonstrations 
of the strong coupling regime in solid-state cQED [3] 
revealed features in the emission spectra for large QD- 
cavity detuning that could not be explained by standard 
Markovian theory [6] . Since then there has been a lively 
debate 0,13:13 ^he origin of the deviations. We demon- 
strate that non-Markovian phonon processes play an im- 
portant role for solid-state cQED. 

Here, using a simple physical model we show that 
photon-assisted dephasing processes are of great impor- 
tance in describing the effect of phonons in a cQED set- 
ting. The underlying physical picture is that the polari- 



ton quasi-particle, formed by dressing the QD with the 
cavity photon, is dephased by phonon processes. We fo- 
cus on the regime of relatively small QD-cavity detunings 
and pulsed excitation conditions where dephasing pro- 
cesses mediated by longitudinal acoustic (LA) phonons 
are expected to be important, and investigate the con- 
sequences on the dynamical properties of the cQED 
system. Pulsed excitation is required for on-demand 
photon sources emitting indistinguishable single photons 
and entangled-photon pairs [ij, making it an important 
regime to investigate. Our theory takes into account 
memory effects of the phonon bath, which are neglected 
in the usual Markovian Lindblad theory of dephasing pro- 
cesses. The non-Markovian formulation is shown to be 
vital in interpreting recent experiments Q- 

A similar non-Markovian theory has recently been used 
to describe the field dependent dephasing dynamics of 
classically driven two- level systems without a cavity Q. 
By varying the strength of the applied classical field one 
can approach a regime where the QD dynamics takes 
place on a time scale that is near or even below the cor- 
relation time of the phonon reservoir, where the usual 
Markovian Lindblad theory of decay breaks down. In 
cQED the coupling is mediated by a single photon with 
a coupling strength g. State-of-the-art samples [7] have 
coupling strengths up to % = 150 /xeV, translating to a 
characteristic time scale of about 14 ps, which is consid- 
erably longer than the phonon reservoir correlation time 
of typically 3-5 ps. Importantly, even for these realis- 
tic parameters non-Markovian effects are found to play 
an important role, giving rise to non-trivial detuning de- 
pendent dynamics and coupling strength. The QD-LA 
phonon interaction and its effects on spectra and dy- 
namics are well understood in the semiclassical regime 
fol, whereas effects due to quantized light fields have not 
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received much attention. Initial work on the influence of 
phonon dephasing on a cQED system has been reported 
[lO| . Common to these works is that httle physical in- 
sight into the non-Markovian nature of the dephasing 
processes has been given, which is a central theme of this 
manuscript. We note that the physical processes identi- 
fied here are also expected to be relevant to other cQED 
systems coupled to vibrational reservoirs, e.g. quantum 
wells, organic molecules, ni trog en vacancy centers in di- 
amond, and colloidal QDs ll], implying that the ideas 
presented here may have wide applications. 

We describe the effect of LA phonons on the QD-cavity 
system using the Jaynes-Cummings model with the ad- 
dition of the electron-LA phonon interaction 12, [l3| . 
The QD has an excited and a ground state of energies 

hio, and hwg, respectively, while the cavity photon has Phonon induced terms in Eq. ^ are 
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FIG. 1: (color online). Evolution of the excited state popu- 
lation sii{t) for hg = 140 /^eV, Iik = 100 /leV, Q ^ 10"', and 
r = 4K. 



energy fiwcav The QD-cavity system space is spanned 
by the two-level basis {|1) — \e,n — 0), |2) = \g,n— 1)}, 
where n is the cavity occupancy, thus |1) describes the 
excited QD and |2) describes the excited cavity. The 
Hamiltonian has three terms, H = Hs + Hi + -ffph, de- 
scribed below. The Hamiltonian of the QD-cavity system 
is Hs = hAaii+hg{ai2+(T2i), where A = Wc— Wg— Wcav is 
the QD-cavity detuning and cTnm = The electron- 

phonon interaction is Hi = an -^^'"(^Ife + ^fe) j where 
M'' = — is the effective phonon interaction ma- 
trix element and creates a phonon in mode k. The 

free phonon Hamiltonian is iJph — X^fc^fc^I^fc' where 
LOf, — Csk is the phonon dispersion and Cg is the speed of 
sound. 

We apply the timeconvolution-less approach [l5| for 
describing the reduced dynamics of the QD-cavity sys- 
tem. The density matrix of the QD-cavity system, 
is considered to second order in iJj, and assuming that 
the phonon bath remains in a thermal state. The equa- 
tions for Spq[t) = Tr{p(i)crpg} (representing the excited 
QD state population, sii{t), the number of photons in 
the cavity, S22(^), and the photon-assisted polarization, 
si2{t)) are 

dtsii{t) = -ig [si2{t) - sl^{t)] - rsii(i), (la) 

dtS22{t) = ig [Sl2{t) - S*12{^)] - KS22it), (lb) 

dtsi2{t) = iAsi2{t) - ig [sii(t) - S22(i)] 

- l/2(r + K)si2{t) + atSi2(i)|ph. (ic) 

Without the phonon induced terms represented by 
dtsi2{t)\ph these equations are the standard lossy Jaynes- 
Cummings model. The losses have been introduced 
through the Lindblad formalism [15|. These include de- 
cay, described by a rate F, of the excited QD state to 
modes other than the cavity and non-radiative chan- 
nels, and the finite linewidth of the cavity, k — Wcav/Q, 
where Q is the usual quality factor of the cavity. We 
take F = 1 ns""'^ in all simulations to be presented. The 



9tsi2(i)|ph = - [7i2(i) - iApoi] si2(i) 

+ ig<{t)s22it)-ig>{t)sn{t). (2) 

This term introduces two novel effects compared to stan- 
dard Markovian cQED models 0]. Firstly, 712 (t) enters 
as a time-dependent pure dephasing rate. Secondly, the 
functions G^{t) renormalize the bare coupling strength 
g: the effective value of g is changed by the real part 
of G^{t), and an additional decay of the polarization is 
induced by the imaginary part of g<{t). The long-time 
polaron shift Apoi = Im {712(00)} has been subtracted 
from 712 (i) [l3|. Explicitly, 



g^{t) = ih-^ f dt'U*,{t')U2iit')D^{t'), 
Jo 



(3) 



712(0 = I dt'[\Unit')\'D<it') - \U2iit')\'D>it'^ 



(4) 

The phonon bath correlation functions are 

D^it) - ^|M'=|2 [uke^'^^' + {uk + 1) eT*"-"*] , (5) 

k 

with Tifc — l/{exp(hu!k/kBT) — 1). The operator U{t) — 
cxp{—iHst/h) is the time evolution operator for the QD- 
cavity system and is the essential ingredient as it intro- 
duces the photon-dressed QD into the phonon scattering 
terms. The Markovian Lindblad formalism is obtained 
neglecting all memory effects associated with the phonon 
interaction, i.e. D^{t) oc S(t) [15|]. In this limit g^{t) = 
and 712 (i) — constant, with no dependence on photon 
properties. 

We have numerically solved Eqs. ([Ij with the initial 
condition of a single excitation on the QD, sii(O) = 1, 
and all other elements set to zero, modeling an experi- 
ment where a QD is excited in its discrete states with 
a short optical pulse [lB|. Theparameters are chosen 
similar to recent experiments [7[. The computed ex- 
cited state populations, sii(i), are shown in Fig. [T]for 
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FIG. 2: (color online), t — > oo limit of ft7i2(t) (blue) and 
hQ^{t) (green/red) for the parameters hg = 140 fieV and 
T — 4:K. Im[ft7i2(cxD)] has the constant value 32.7 fieV. 
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FIG. 3: (color online). Real part of ^712(00) (solid) and 
hGRioo) (dashed) for Fig = 50 (blue), 100 (red), 150 (green) 
fieV. The imaginary parts are independent of temperature 
and much smaller than the real part and therefore not shown. 



various detuning values. For A = we observe an ex- 
pected strong enhancement of the decay rate and asso- 
ciated Rabi oscillations, indicating the strong coupling 
regime. For non-zero detuning we observe small Rabi 
oscillations in the start of the decay curve, again indica- 
tive of the strong coupling regime. Interestingly, our the- 
ory predicts a shorter lifetime when the cavity is tuned 
below the QD resonance than when it is tuned above, 
in stark contrast with standard Markovian theory [l7| . 
The asymmetry is particularly strong for the A/k = ±10 
cases and has recently been observed in experiments 0. 
We attribute this pronounced reduction in lifetime to a 
phonon-assisted Purcell effect, where the QD may cou- 
ple via phonon emission to the cavity when the cavity 
is spectrally below it. In the opposite case, where the 
cavity is spectrally above the QD resonance, phonon ab- 
sorption is needed for the QD to become resonant with 
the cavity. This process is suppressed at low temper- 
atures and the phonon mediated coupling between QD 
and cavity is lost. The asymmetry becomes less pro- 
nounced when the detuning is larger than the energy the 
electron can lose through phonon emission, as seen in the 
A/k — ±30 cases. The importance of using the photon- 
assisted electron-phonon interaction, (i.e., the operator 
U{t) in Eqs. ([3]) and (|4])), should be emphasized. Set- 
ting g = only in U{t) we obtain results in quantitative 
agreement with the A < results in Fig. ([ij, where 
the phonons do not significantly influence the dynamics, 
for all values and signs of A. Formally this is easily un- 
derstood, as in this case G^{t) = and 712(0 loses its 
dependence on g and A, and thus the asymmetry is lost 
in the phonon induced dephasing. One can interpret this 
approximation as only allowing the phonons to interact 
with the bare electron and not the electron-photon quasi- 
particle, the polariton, that is actually present in the 
system. This illustrates the importance of accounting for 
the polaritonic quasi-particle nature of the strongly cou- 
pled QD-cavity and its non-Markovian interaction with 
the phonon reservoir. 

Figure [2] shows the long-time limit [3l of the rates. 



Eqs. dni) and A strongly asymmetric Q^{oo) is ap- 
parent versus detuning. Also, ty^(oo) and ^?^(oo) attain 
very different values for fixed detuning. As the pure de- 
phasing rate 712(00) remains symmetric with respect to 
detuning, the asymmetry in Fig. [T]is caused by the renor- 
malization rates Q<(oo). It is remarkable that despite of 
their relative weakness, |^^(oo)|inax/3 < 7 %, such large 
dynamical effects may occur. The physical origin of this 
asymmetry can be traced back to the phonon correla- 
tion functions given in Eq. ([5]). Here it is seen that 
phonon absorption processes are suppressed at low tem- 
peratures as they are proportional to the occupation fac- 
tor rtfc, whereas phonon emission processes continue to be 
possible due to the presence of the phonon vacuum field. 
In linear QD absorption spectra virtual phonon emission 
results in asymmetric sidebands centered around an in- 
finitely sharp zero phonon line (ZPL) [19]. However, for 
the present physical system it is essential that the cav- 
ity field is treated non-linearly. This results in novel ef- 
fects such as a finite width of the ZPL and the coupling 
strength renormalization rates G^{t). 

We next examine temperature effects at zero detuning. 
For this case: g^{t, A = 0) = C?ij(t) T«5/(0> where gR{t) 
and Giit) are real functions. Figure [3] shows the effect of 
temperature on the phonon induced rates Q^{oo) and 
712(00), within the low temperature regime typically ex- 
plored in cQED experiments. As expected, the pure de- 
phasing rate increases with temperature. We also observe 
an increase as a function of the bare coupling strength 
9- S_r(oo) also increases in magnitude with increasing 
temperature, but has a negative value. This leads to a 
lowering of the effective coupling strength entering the 
equation for si2{t) as temperature is increased. This can 
be realized by inserting Eqs. ^ and ^^(00, A = 0) into 
Eq. He]) yielding geff(oo) = 5 - |^fl(oo)| < g. The ef- 
fective coupling strength thus depends significantly on 
temperature. This mechanism is also relevant for inter- 
preting cQED absorption spectra [S]- We expect this 
dependence to have a detrimental effect on the possibility 
of reaching the strong coupling regime. To quantify this 
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FIG. 4: (color online). Parameter space showing the presence 
of strong or weak coupling for the full model (solid) and for 
G'^it) = (dashed), with Hk = 75 (red), 125 (blue), 175 
(green), 225 (black), 275 (magenta) /xeV. 

prediction we have investigated the transition between 
the weak and strong coupling regime, while varying the 
most important parameters in the model, namely g, k, 
and T. We define the weak coupling regime as the situa- 
tion where the initially populated excited state of the QD 
decays monotonically toward zero. Figure H] shows the re- 
sults both for the full model and for comparison the case 
where we have artificially put G<{t) = 0. This is done to 
emphasize the effect of the temperature induced renor- 
malization of 5, by allowing only the pure dephasing rate 
to be temperature dependent, as is common practice in 
phenomenological cQED models. As expected, we gener- 
ally observe the presence of strong coupling in the system 
for large g and low T, with the parameter space of strong 
coupling becoming extended as we increase the quality of 
the cavity. Comparing the results of the full model with 
those where G<{t) = we notice a strong effect of the 
renormalization of the bare coupling constant g. The 
parameter space where strong coupling is obtained is sig- 
nificantly decreased when including the renormalization 
of g. This result is relevant in the interpretation of ex- 
perimental data, as state-of-the-art cQED models @] do 
not include the renormalization effects contained in the 
functions G<{t). These effects are of significant impor- 
tance and therefore cQED models neglecting them can 
lead to misinterpretation of experimental data. 

In conclusion, we have illustrated the importance of 
applying a dressed-state picture of the polaritonic QD- 
cavity system when modeling cQED systems interacting 
with LA phonons. We have shown its relation to recent 
experiments, explaining the observed asymmetry in life- 
times with respect to the QD-cavity detuning. The asym- 
metry can only be understood by treating the QD-cavity 
system as a polaritonic quasi-particle in the phonon in- 
duced scattering terms. Furthermore we have investi- 
gated a phonon induced lowering of the effective coupling 
strength with increasing temperature, which was found 
to change the criterion for strong coupling significantly. 



While completing this manuscript we became aware of 
a recent preprint with related work |20| . however with 
focus on modeling experimental emission spectra as op- 
posed to interpretation of the mechanism responsible for 
the dynamics, which is presented here. 
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